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Abstract 

A detailed chemical dynamic model is presented for a moist Martian 
atmosphere. Recombination of CO 2 is catalyzed by trace amount of H 2 O. 

The abundances of CO and Og should vary in response to changes in 
atmospheric H 2 O and atmospheric mixing. 
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Carbon dioxide is the major gaseous component in the atmospheres 
of both Mars and Venus and it is difficult to understand its apparent 
stability,. The gas is readily dissociated by sunlight at wavelengths 

O 

less than 2000A. However, recombination of C0 2 by the elementary 
reaction 

CO + 0 + M - C0 2 + M (1) 

proceeds at a rate which is negligibly slow (1) compared with the competing 
reaction 

0 + 0 + M - 0 2 + M (2) 

It would appear therefore that the atmospheres of Mars and Venus should 
contain large amounts of 0 2 and CO, and yet both atmospheres are remarkably 
deficient in dissociation products such as CO, 0 2 , 0^ and 0. 

We shall argue here that recombination occurs predominantly between 
CO and 0 and proceeds catalytically by the reaction sequence 

H + 0 2 ,+ M H0 2 + M (3) 

0 + H0 2 "* OH + 0 2 (4) 

CO + OH =* C0 2 + H . (5) 

Photolysis of H2O provides the source of hydrogen radicals (2). This 
sequence readily accounts for removal of oxygen atoms below about 25 km. 

At higher altitudes the H0 2 concentration is vanishingly small, a con- 
sequence of the strong dependence on altitude of the three body reaction (3). 
The alternate schone 


H + O3 - OH + 0 2 


(6) 
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followed by (5) plays some role at higher altitudes where ozone is produced 
mainly by 

0 + 0 2 + M “• 0^ + M (T) 

and removed by photolysis in the Hartley continuum 

hv + 0g “» 0 + 0 2 . (8) 

A viable photochemical model must also account for a chemical 
balance in 0 2 . Molecular oxygen is formed by (2) and by 

0 + OH -» 0 2 + H (9) 

Both reactions provide large potential sources of 0 2 above 20 km. Molecular 
oxygen diffuses downward and is removed by photolysis in the Herzberg 
continuum 

hv + 0 2 -* 0 + 0 , (10) 

followed by reactions ( 3 ) - (5)> with some additional removal associated 
with the reaction sequence ( 3 ) 

HO 2 + ho 2 - h 2 o 2 + o 2 (11) 

hv + H 2 0 2 -* 20H (12) 

followed by reaction ( 3 ). However, the bulk of the recombination is due 
to reactions ( 3 ) - ( 5 ). This is a consequence of rapid downward transport 
of oxygen atoms which suppresses the formation of 

An alternate scheme (2,4) for recombination of CO and 0 2 involves 
reaction ( 3 ) followed by 


CO + H0 2 - C0 2 + OH 


( 13 ) 
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and reaction (5) . We shall argue that reaction (13) is too slow to play 
a role for Mars. 

The importance of HgO in the Martian atmosphere was emphasized 
earlier in another context (5 ) . Recombination of 0 2 + in the Martian 
exosphere leads to production of energetic oxygen atoms and a significant 
number of 0 atoms escapes to interplanetary space. The estimated escape 
rate for 0 is approximately equal to half the observed escape rate for H. 
The H atoms are supplied to the upper atmosphere by upward diffusion of 
Hg formed as a by-product of H 2 0 chemistry in the near surface region. 

The primary source of H 2 is the reaction 

H + H0 2 >- H 2 + 0 2 (14) 

and the net production of H 2 is determined by the magnitude of the 0 atom 

escape rate. It was argued that the relative concentrations of CO and 0 2 

in the lower atmosphere would be moderated as a result of oxygen escape 

so as to supply the required upward flux of H 2 . Escape of oxygen acts 

therefore to determine the magnitude of the quantity 

00 

1 = £ k l4 [H][H0 2 ] dz (15) 

where [h] and [HOg] are number densities of H and H0 2 at height z and 
k^ 4 is the rate constant for reaction (14). Using the analyses by Hunten 
and McElroy (2) and McElroy (5) we estimate I = 3*5 x 10® cm"^ sec"- 1 - with 
a probable uncertainty of about a factor of 2. An acceptable chemical 
model for Mars is constrained therefore to provide the appropriate value 
for 1.(6) 

Detailed numerical studies of Martian photochemistry are simplified 
by a number of factors. First, O 2 and CO may be assumed to be in diffusive 
equilibrium over an extensive altitude regime. Hie relevant photochemical 
time constants are long compared to any conceivable transport times. We 
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may therefore assume that the abundances of CO and Og are known and adopt 
the observed mixing ratios of 8 x 10”^ (7) and 1.3 x 10“3 (i3) respectively . 
Secondly, the concentration of odd hydrogen, [h] + [HO 2] + [OH], should 
also be well mixed in view of the relatively long time constants for 
radical removal. The principal paths for radical loss are (14) and 

OH + H0 2 =* HgO + 0 2 (16) 

The dominant radical below 25 km is H0 2 . Atomic hydrogen dominates at 
higher altitudes and the mixing ratio of odd hydrogen, for given Og and 
CO, is specified by the conservation condition (L5). With the model 
employed here we find a mixing ratio for odd hydrogen equal to 5-0 x 10“ 10 . 

Hie detailed chemical model is summarized in Table 1. The important 
rate constants are moderately well known, with two exceptions. Reaction (11) 
is currently in some dispute. Clark (l), noted that extrapolation of high 
temperature data obtained by Baldwin et al. (9) implied a rate constant 
at Martian temperatures of order 10“30 cm3 sec"^. Davis (10) determined 
an upper limit for the reaction at room temperature equal to 10“-^ cm3 
see”' 1 '. On the other hand Westenberg and DeHaas (11) concluded that the 
rate constant at room temperature was fast, equal to approximately 0.06 
times the rate constant for reaction (15). The choice of rate constant 
in Table 1, k^ < 10 =1 ^ cm3 sec” i , reflects a careful study of the 
available data. We note for example that Westenberg and De Haas (11) 
did not allow for the fast reaction (16) as a loss mechanism for OH. On 
the basis of our current understanding of Mars it would be difficult to 
accept a value of k^g much larger than the upper limit adopted here (12). 

Same uncertainty is associated also with k^. Kauftaan (13) con- 

<=12 3 -1 

eluded that k-^ should exceed 3 x 10 cm- 5 sec , 


reflects a recent 
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analysis of the available data by McConnell (14) . A smaller value for 
ki^ would imply somewhat larger mixing ratios for odd hydrogen and lower 
values for the eddy mixing coefficient near 30 km but should not signifi- 
cantly alter the general characteristics of the photochemical model. 

In particular, recombination of 0 and CO would continue to occur primarily 
by reactions (3) - (5), although photolysis of H 2 O 2 would be relatively 
more important as a sink for O 2 . 

Concentrations of odd oxygen were obtained by numerical solution 

of the coupled diffusion - continuity equations. Odd oxygen is produced 

by photolysis of C0 2 and 0 ^. We used recent data by Widing et al. (15) 

to compute dissociation rates. Their fluxes are somewhat lower than 

values employed in earlier calculations (l6) and photolysis rates are 

similarly reduced. The net production of 0 in the present model is equal 

to 1.8 x 1012 cm -2 sec” 1 . Loss of odd oxygen is primarily by (5), with 

some contribution from (9) and (2). Results for the principal constituents 

are shown in Figure 1. This model corresponds to a vertical eddy diffusion 

a o 

coefficient of 1.5 x 10 cm sec and is most sensitive to the diffusion 
coefficient in the vicinity of 30 km. The magnitude of the diffusion 
coefficient derived here is similar to values discussed by Gierasch and 
Goody (17) in their study of energy transfer in the Martian troposphere. 
According to calculations by McElroy and McConnell (16) somewhat larger 
eddy coefficients, ~ 10^ cm 2 sec”' 1 ', are required at higher elevations 
(above 60 km) in order to account for upper atmospheric data on 0 and CO. 

The relative importance of various source and sink terms is 
illustrated in Figure 2. The quantities plotted here are rates of various 
reactions, integrated from the Martian surface to height z (km). The 
integrated rate of CO 2 photolysis is also shown. The difference between 
this rate and the total recombination rate reflects the contribution of 
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flow to the local atomic oxygen production rate. Evidently CO-O recom- 
bination occurs mainly below 25 km, and the reaction sequence (3) - (5) 
is the dominant path for recombination. Formation of Og occurs in a 
restricted altitude regime, between 25 km and 30 km, and reaction (9) 
dominates. The rate of reaction (l6) obtained from Figure 1 implies an 
average rate for photolysis of 1^0 equal to 2.7 x 109 cm" 2 sec"^, in 
agreement with earlier estimates (2) of the mean photolysis rate based 
on observed H2O concentrations (18) . The model predicts an O3 abundance 
of 1.4 x 10"^ cm atm., in satisfactory agreement with limits set by the 
ultraviolet experiments on Mariners 6, 7, and 9 (19) • Ozone Is normally 
a minor source of ultraviolet opacity in the Martian atmosphere although 
occasionally strong absorption is detected. High ozone concentrations 
are apparently correlated with unusually cold atmospheric conditions, (19) > 
and may be a natural consequence of the chemical model discussed here. 

One would expect lower concentrations of H 2 0 and consequently OH, H0 2 and 
H in colder regions of the Martian atmosphere. Wet chemistry is less 
effective and odd oxygen concentrations should be consequently higher. 

This matter will be discussed in detail elsewhere. 

The response of the Martian atmosphere to a change in the concen- 
tration of H^O is of considerable interest. For modest changes, the 
atmosphere readily adjusts to a new equilibrium state with a different 
mixing ratio of CO relative to 0 2 . To first order, the concentration 
of OH is proportional to the concentration of H2O and the concentrations 
of 0, H, and H0 2 in the lower atmosphere remain constant. Only CO 
changes, and the change is inversely proportional to the change in H 2 0. 

To a higher level of approximation, however, one must take account of 
consequent changes in the 0 2 production and one finds necessarily some 
decrease (or increase) in 0 2 associated with a decrease (or increase) 
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in HgO. By numerical experimentation we found that if the HgO concentra- 
tion were arbitrarily reduced by a factor of 2, the CO concentration would 
increase by a similar factor and 0 2 would decrease by about 50 $, if the 
eddy mixing coefficient remained constant. 

It is of interest to consider also the atmospheric response to 
a change in the eddy mixing coefficient near 30 A reduction (or 

increase) in the eddy mixing coefficient at 30 km leads to an increase 
(or decrease) in 0^ formation. The chemical equilibrium in the lower 
atmosphere is altered with a consequent increase (or decrease) in the 
abundance of CO and 0 2 . The mixing ratio of CO relative to 0 2 remains 
relatively constant. The suggestion of a time variable 0 2 concentration 
in the late stages of the recent planet=wide dust storm (8) may reflect 
in part a change in the dynamic state of the atmosphere, in part a change 
in atmospheric 1^0. It is difficult however, to draw more quantitative 
conclusions. A complete analysis should allow for heterogeneous chemistry 
on atmospheric dust and there are reasons to suspect that Martian dust 
may play some role in the chemistry of atmospheric oxygen (20). 

Finally, we note the central importance in the present analysis 
of the conservation condition afforded by equation (15). The relative 
abundance of odd hydrogen, and the mixing ratio of odd hydrogen, are 
controlled by this relation in the present scheme. An increase (or decrease) 
in 0 escape, associated for example with increased (or decreased) solar 
activity, will result in an increase (or decrease) in atmospheric CO. The 
concentrations of OH in the lower atmosphere will change. For fixed HgO, 

HO 2 will increase (or decrease) corresponding to an increase (or decrease) 
in 0 escape. Ihe Hg production and H escape will alter accordingly. The 
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present model provides plausible justification for the observed relation 
between H and 0 escape and lends confidence to conclusions derived 
earlier (5) regarding the evolutionary history of Martian H^O (23). 

Michael B. McElroy 
Thomas M. Donahue 

Center for Earth and Planetary 
Physics, Harvard University 
Cambridge, Massachusetts 02138 
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Table 1 

Relevant reactions with their rate constants. For 2 body reactions, 
units are cm3 sec”- 1 - and for 3 body reactions, units are cm^ sec"- 1 -. 


1 

co + o + co 2 - co 2 + co 2 

k x = 2 x 10" 37 

2 

0 + o + co 2 o 2 + co 2 

kg = 3 x 10 " 33 ( T / 3 oo )“ 2 ‘ 9 

3 

H + 0 2 + C0 2 -» H0 2 + C0 2 

k^ = 2 x 10” 31 (t/273 )“ 1 * 3 

k 

0 + H0 2 -* OH + 0 2 

k^ = 7 x 10" 11 

5 

CO + OH C0 2 + H 

k^ = 9 x 10“ 13 exp(~500/T) 

6 

H + O 3 OH + 0 2 

k^ = 2.6 x 10" 11 

7 

0 + o 2 + co 2 ““03 + co 2 

k 7 = 1.4 x 10 " 33 ( T / 3 oo )" 2 

8 

0 ^ + hv -* 0 2 + 0 

Jg = 4.2 x 10 " 3 see " 1 

9 

0 + OH -* 0 2 + H 

k 9 = 5 x 10" 11 

10 

0 2 + hv -* 0 + 0 

J 10 = 5.8 x 10" 10 sec” 1 

11 

CO + H0 2 - C0 2 + OH 

ku = < IO " 16 

12 

ho 2 + ho 2 -» h 2 o 2 + o 2 

k^ = 9.5 x 10" 12 

13 

H 2 0 2 + hv - OH + OH 

J 13 = 5.2 x IO" 9 sec" 1 

14 

H + H0 2 -• Hg + 0 2 

k l4 = 1 x IO" 11 

16 

OH + HOg - HgO + 0 2 

k l6 = 2 x IO" 10 


General References; Clark (l), Kauftnan (13) 

Rate Constants: k x (22, 23, 24), 1^(25), ^(21), 1^(26), k 13 (27), k l6 (2l). 
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Fig, 1 


Fig. 2 


Figures 


Concentration of principal constituents in the Martian 
Atmosphere. Surface temperature is 220°. Eddy diffusion 
coefficient is 1.5 x 10° cm see” . The curve for odd 
hydrogen, shown as a broken line between the portions 
labelled HOg and H, is plotted at 5 x 1CT 10 times the 
CCS£ density. The 0g and CO densities are respectively 
1.3 x 10“3 and 8 x 10"^ times the COg densities. 

Integrated reaction rates important in COg and Og 
formation and loss. The curve labelled P(o) is the 
integrated photolysis rate for COg. The contribution 
from 0 + 0 + M is only 2.3 x 10^ cm 2 sec =1 and not 


shown 






Figure 2 




